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Abstract. The paper develops a general theory for finite rubber viscoelasticity, and 
specifies it in the form, convenient for solving problems important for rubber, tire and 
space industries. Based on the quasi-linear approach of non-equilibrium 
thermodynamics, a general nonlinear theory has been developed for arbitrary non- 
isothermal deformations of viscoelastic solids. In this theory, the constitutive equations 
are presented as the sum of known equilibrium (rubber elastic) and non-equilibrium 
(liquid polymer viscoelastic) terms. These equations are then simplified using several 
modeling arguments. Stability constraints for the proposed constitutive equations are 
also discussed. It is shown that only strong ellipticity criteria are applicable for assessing 
stability of the equations governing viscoelastic solids. 


1 INTRODUCTION 


Within the general rational mechanics 1 , two types of continuum approaches to nonlinear 
constitutive equations (CE’s) for viscoelastic solids have been suggested for applications. 
One is of the differential type (e.g. see 2 ’ 3 ), well developed for viscoelastic liquids on the 
basis of non-equilibrium thermodynamics. The other approach is the K-BKZ ’ single- 
integral CE’s. Both of them might be considered as complimentary in application to 
viscoelastic solids, and more precisely, to the thermo-mechanical and relaxation 
properties of cross-linked rubbers. 

The effects of relaxation phenomena, which do not exist in pure elastic solids, 
enhance complexity of problems when polymers are considered. Hence, viscoelastic 
polymer solids demonstrate many nonlinear rheological effects, such as nonlinear creep 
and relaxation, hysteresis etc., which are not observed in other solids. Because of this, 
even geometrically simple problems that are analytically solvable for pure elastic solids 
have to be treated numerically in the viscoelastic case. 

There are practical needs for developing a thermodynamically based approach to 
solid nonlinear viscoelasticity, which would be valid over a wide range of strains, strain 
rates (frequencies) and temperatures. It seems that this approach should combine (i) 
nonlinear multi-mode CE’s of differential type generated by a part of relaxation spectrum 
with high relaxation times, and (ii) nonlinear single integral CE generated by that part of 
relaxation spectrum with low relaxation times. We suppose that both parts of the 
combined CE are almost independent. Part (i) governs the relatively slow motions of a 
continuum undergoing very large strains and obeying the WLF time-temperature 
superposition principle at relatively high temperatures. Part (ii) governs the fast (up to 
very fast) motions of continuum undergoing relatively small (even very small) strains and 
obeying a different time-temperature scaling that is valid for the glassy state. Thus the 
approach takes into account thermo-rheological complexity along with the high non- 
linearity of elastic and viscoelastic deformations, including compressibility effects. To 
make the approach reliable, the constraints of strong hyperbolicity and thermodynamic 
stability are imposed on the developed CE’s. In this paper, only the first part of the above 
proposal, the part (i) of the CE’s, is demonstrated. The second part (ii) will be displayed 
elsewhere. 

In models of a differential type for viscoelastic solid (cross-linked) rubbers, the 
typical nonlinear pure elastic term should be included along with differential modes for 
different and well-separated relaxation times. The simple assumption is that this term is 
additive to the viscoelastic ones. In typical development 6,7 of these models, only a 
nonlinear extension of the “standard” (three-parametric) linear viscoelastic solid model 
has been used. Also, the treatment of the dissipative term in the evolution equation for 
nonlinear Maxwell element in paper 6 is not convincing, and the whole approach taken in 
the development of the evolution equation for this element in paper 7 is questionable. An 
additional important question to ask of paper 7 , where both the pure gum elastomers and 
carbon black filled compounds were studied, is its treatment of Mullins hysteresis for 
carbon black filled rubbers as a pure relaxation effect. It is known, however, that this 
effect, being structurally reversible (i.e. thixotropic), takes several months at room 
temperature for a complete recovery 8 . The differential modeling has another importance: 



one can establish the transformation of the CE from liquid to solid-like behavior via 
vulcanization using chemorheological modeling 9 . Surprisingly, this type of modeling was 
done first for filled elastomer compounds 8 . 

Other nonlinear viscoelastic models of differential type have also been recently 
proposed in the literature, being based on pure continuum mechanics phenomenology. 
For example, several models have been developed (see paper and references there) that 
include nonlinear elements for modeling plasticity effects with a von-Mises type yield 
criterion. Also, in paper 11 the authors attempted to make a nonlinear extension of the 
linear 3-parametric “standard” model with two equivalent constructions of springs and 
dashpot. The two nonlinear models have been presented in an awkward form, based on 
the separation ‘'the total viscoelastic strain into a product of elastic and inelastic 
components, with the use of the Second Law. Another approach , close in spirit to that 
elaborated here, has been recently introduced. There the author analyzed the viscoelastic 
solids in terms of a separated set of nonlinear Maxwell modes (with hidden tensor 
parameters), added to the equilibrium elastic mode. Although the equilibrium mode has 
been treated in a general and correct way, the evolution equations, as written, violate 
material objectivity, and the dissipation term was taken to be similar to the stress in the 
Maxwell mode (i.e. without strain induced anisotropy). The same defects in analyzing 
viscoelastic evolution equation in models for finitely deforming viscoelastic solids 
(extended from the three-parametric “standard” solid viscoelastic model) are observed in 
papers 13 ' 15 . Stability constraints were not discussed in papers " . 

The first step in a practical development of nonlinear models for viscoelastic solids 
using a multi-modal differential approach is elaboration of reliable methods for 
determining linear relaxation spectra. Since modern molecular models (e.g. see ' ) are 
still not reliable, the relaxation spectrum must be determined from experimental data. 
Because this problem is ill-posed, some mathematical methods have been elaborated. The 
Pade-Laplace method 19 ’ 20 , developed for viscoelastic liquids, is the easiest and most 
reliable tool for effective discretization of relaxation spectra by a Prony series. When 
applying the Pade-Laplace method to a viscoelastic solid (e.g. for cross-linked elastomer) 
the true elastic component should be found first, via a preliminary experiment, and then 
removed from the Pade-Laplace procedure. 2 3 

The structure of the paper is as follows. Adopting the thermodynamic approach ’ 
developed for nonlinear liquid-like viscoelasticity, we initially develop a general non- 
linear theory for elastic and viscoelastic nonlinear behavior of solids, with a single 
nonlinear relaxation mode. We simplify the developed theory and then extend this 
simplified format for the multi-mode case, along with a simplified fractioning of 
compressible and incompressible effects in nonlinear elasticity and viscoelasticity. We 
finish up by discussing the stability constraints imposed on the general theory, as well as 
on the simplified and extended models. 

2 STATE VARIABLES AND FREE ENERGY 

Here we will employ the general approach of irreversible thermodynamics discussed in 
various books (e.g. see 21 ). We assume that in the simple case under study, the state 
variables for viscoelastic solid are: the temperature T, a measure of total strain, say, the 
Finger tensor B , and a set of hidden variables that are taken to be symmetric, positive 



definite, second-rank, non-dimensional tensors b . Each tensor b is treated as the 

—k —k 

Finger elastic (recoverable) strain, with the capability of being independently measured. 
It is expected that a tensor b is associated with the nonlinear k tb relaxation mechanism 

that, can hopefully can properly describe the (highly) nonlinear relaxation phenomena 
observed in cross-linked rubbers. In the following, we employ the Eulerian formulation 
for presenting our constitutive equations and, without loss of generality, use a Cartesian 
coordinate system. 

To start off, we consider a single relaxation mechanism, whose elastic energy 
associated with a non-steady deformation is described by the tensor b . Being mostly 

interested in studies of isothermal deformation of isotropic viscoelastic solids, we 
introduce as a proper thermodynamic potential, the Helmholtz’ free energy density F per 
mass unit, which depends on temperature T and the basic invariants of tensors 5 and b. 

The crucial simplifying assumption underlying all of today approaches is that the free 
energy might be represented as a sum of two contributions, one from the equilibrium 
processes (with tensor# as its parameter), and another from non-equilibrium processes 

(with tensor b as its parameter): 

F~F\T,I?,I b ,I b ) + F 1 (T-I v I 2 ,I 2 ). (1) 

Here I B and l - are respectively the basic invariants of tensors # and b , defined as: 


I B =trB, I B = \/2(I 2 -trB 2 ), / B =det#; 

/, = trb, 1 2 = U 2(1 ; 1 - trb 2 ), / 3 = det £ 

Then the “thermodynamic” stress tensor cr r associated with the free energy F, is: 

£ r== 2 0 + £i’ 

wherein 

£ o = 2pB dF° /dB = 2p[F° B + F° (I? B-B 2 ) + F 3 °I 3 b S] ( F°=dF°/dI B ) 

(4) 

a =2pbdF'/db = 2p[F 1 , b + F 2 l (/ 1 b-b 2 ) + F 3 'l 3 £]. (Fj=dF/dl\) 

In eqs.(3) and (4), p is the density, § is the unit tensor, £ o and q _ ( are pure elastic 

(equilibrium) and viscoelastic (quasi-equilibrium) stresses, respectively. The only 
physical reason to use the suggested form in the second equation in (4) is the explicit 
assumption 2,3 that polymeric liquids and solids always possess an “instantaneous” elastic 
limit, a quasi-equilibrium situation achieved on very fast (instantaneous) deformations, 



where the temporary entanglements in macromolecules act like additional cross-links. 
Note that no incompressibility has been assumed here. 

The well-known kinematical relation for the Finger total strain 5 , 

v o 

B = B-Be-eB=0, (5) 

v o 

will be used in the following, where 5 and 5 are upper convected and corotational 
(Jaumann) tensor time derivatives, and eis the strain rate (or rate of deformation) tensor. 

The existence of a like kinematical (or kinetic ) equation that governs the evolution of 
tensor b , the main goal of the following constitutive theory, has yet to be established. 

3 ENTROPY PRODUCTION AND NONLINEAR NON-ISOTHERM AL 
RELAXATION PHENOMENA 

2 3 

The dissipative effects have been thoroughly analyzed for viscoelastic liquids^' . In 
viscoelastic solids, these effects are associated with their time dependent deformations. 
The Clausius-Duham expression for the entropy production P s has the form: 

TP S = —q ■ VT + tr(a ■ e) - pdF/dt\ T (> 0) . (6) 

Using (4) and (5), we can rewrite eq.(6) as follows: 

. o 

TP s .=-q VT + tr[(a-a o )-e]~ tr (<7 ] -b 1 / 2 b) . ( 7 ) 

The next step is introducing the kinetic relation for tensor b , similar to the kinematical 
equation (5). This has been intensively discussed in Ref. 3 and resulted in the definition of 
a “thermodynamic flux”, the tensor e , such that 

b—be^—e^ e b = 0. ( 8 ) 

Note that in the equilibrium limit, b^>B, —3 ► £ , CT f — > ct q , and eq.(8) coincides with 

eq.(5). Therefore it is convenient to introduce the pure non-equilibrium (dissipative) 
quantities, the stress o_ and strain rate that vanish in the equilibrium: 

= (J-a T - <7 -<J Q -£ J ; (9) 

Substituting the second relation in (9) into eq.(8) reduces the kinetic equation to the form: 



( 10 ) 


v 

b+ b ■ e + e p ■ b = 0 . 

Here the dissipative strain rate tensor e^has yet to be established, and is expressed as a 

function of temperature and the tensors 7? and b . 

Substituting also eqs.(9) into (10) into (7) finally yields: 

TP S - —q - VT + tr{a^ ■ e) + tr( £ ( -e) (>0). (11) 

The entropy production in eq.(43) is now presented as a typical bilinear form 
where) X k }={ VT,a ,q_ x }are the independent thermodynamic forces of different tensor 

dimensionalities, and conjugated to them are the quantities, { Y k }={ q,g,g p }, which are 

independent thermodynamic fluxes. Three independent dissipative sources are now 
clearly seen in equation (11). They are due to the non-isothermality (the first term in 
eq.(ll)), and two contributions in mechanical dissipation in the system (the second and 
third terms in eq.(ll)): (i) the mechanical dissipation produced by work of the 
irreversible stress <7 on the total strain rate g, and (ii) the mechanical dissipation 

produced by work of the viscoelastic stress on the irreversible strain rate e . 

Using the arguments of tensor dimensionality and Onsager symmetry of kinetic 
coefficients, we can write the CE’s in the general quasi-linear form: 

a p =M, :e-M 2 :£,; e p = M 2 : £ + M 3 : a ; M k = M k (T,B,b ) ; (12) 

£ = -£'VT; £ = £( 7 '’!’^) (13) 

Here the kinetic coefficients M k are tensors with rank of four, symmetric within the first 
and second pairs of indices, and symmetric by transposition between the first and second 
pairs of indices. Tensor M, has a dimensionality of viscosity, tensor M 2 is non- 
dimensional, and tensor M 3 has a dimensionality of reciprocal viscosity. The scalar 
products in eq.(12) mean, (M, :g), y =M' jkl e lk ,e tc. The heat conductivity £is presented 

in eq.(13) as a second order symmetric tensor. The dependences of the kinetic tensors in 
eqs.(12) and (13) on the strain tensors B and b demonstrate the effect of strain induced 

anisotropy. Substituting eqs.(12) and (13) into eq.(l 1) presents the entropy production as 
the positively definite quadratic form: 

TP s = tr(K ■ VTVT) + 77? (M, : ee) + TR(M 3 :£,£,) (> 0) . (14) 

Due to eq.(14) the second rank tensor rand the two rank of four kinetic tensors M,and 
M 3 are positively definite. Eq.(14) imposes no thermodynamic constraint on the kinetic 



coefficient M 2 from eq.(12). The complicated general structure of the rank-four kinetic 
tensors was explicitly exposed in Ref. 2 ~. The general structure of the heat conductivity 
tensor k(T,B,c), also quite complicated, can be readily established using the 

representation theorems 1 (p.35). Equations (10), (12) and (13) present the closed set of 
CE’s for non-isothermal, finite, solid, compressible viscoelasticity. 

Introducing the heat capacity c d under constant tensors 5 and b as c d — dU / dT gg , 
reduces the balance of internal energy U to the common heat equation: 

pc d dT I dt = V •£• VT + tr(&_ e) . 05) 

In the case of rubber (entropic) elasticity, when U -U(T), the heat capacity^ is the 
function of only temperature. 

The above general phenomenological relations include in the consideration the 
Kelvin-Voight stress that denies the commonly known “instantaneous” elastic response 
that has been well documented in experiments. Therefore in the following, the only 
admissible type of solid viscoelasticity, the Maxwell type with instantaneous response, 
will be analyzed. This considerably simplifies the above CE’s to the form: 

£ = £ r =£o + £, (M, =0, M 2 = 0) ; e p =M (T,B,b):a (M = M 3 ) . (16) 

Here the expressions for £ q and a are given in eq.(4). 

4 COMPRESSIBILITY EFFECTS 

Although the above general theory includes the compressibility effects, it is convenient to 
expose them explicitly using the density as a new variable. We analyze initially the 
compressibility effects in single-mode viscoelastic approach with instantaneous elasticity. 

In the equilibrium (pure elastic) case, the free energy is represented in the form: 
F = F c {T'jf ,It ,1*) , the stress tensor is given by the second formula in (4), with the 
kinematical relation (5). Additionally, there are two well-known equations for mass 
conservation: 

/ 3 b =(p 0 //?) 2 , d,p + V-(pu) = 0 or d\npl dt = -tre. (17) 

Here u is the velocity vector. Note that the second relation in eq.(17) is easily derived 
from the first one and the kinematical equation (5). Using the first relation in eq.(17) we 
now introduce instead of tensor B , its “incompressible’ part B and density p , 

B = (I^y' n B = (p/p 0 ) in B, 

with corresponding invariants 7 y e of tensor B as: 


(18) 


lf=trB = (plp Q ) 2n lf, i°=trB' = (p/p 0 ) in I* , if = 1. (19) 


Using eq.(19) the free energy F ° , can be expressed via the density p and new invariants 
if and If as: F° ( T , if , if ,lf) = F°(T, p, if Jf ), with the representation of the stress 
tensor as: 


£° = -P 2 F?£ + 2 p[Ff (| - 1 / 3 If S) - Ff ( J~‘ - 1 / 3 if S)] . (20) 

Eq.(20) has been derived by Truesdell 23 . The first term there represents the equilibrium 
thermodynamic pressure, and the second, deviatoric term, demonstrates the effects of 
isochoric (“shearing”) deformations. Finally, using eqs.(17) and (18) reduces eq.(5) to the 
“incompressible” form: 

B = B-Be-e-B=0 ; e^e-l/35tre = e-l/3Sp~'dp/dt. (21) 

The above derivation is repeated in the non-equilibrium, viscoelastic case, starting with 
the introduction of a new irreversible (“density-like”) variable p l as: p' / p Q =/ 3 l/2 , and 
leading to the similar formulae, 

b = (I,y U3 b^(p l /p 0 ) 2n b, (18a) 

/, = trb = (p‘ / p 0 ) 2n I { , i 2 = trb = (p' //? 0 ) 4/3 / 2 , / 3 = 1 . (19a) 

Using eqs.(18a), (19a) reduces the evolution equation (10) to the following: 

o 

b-b e-e b + b e +e b = 0; d\np'/dt-tre = Jin p / dt = -tre . (22) 


Here e and e are the deviators of the tensors e and e ^ , respectively. 

Finally, using for the non-equilibrium part of free energy F ] the presentation 


F\T,I X ,/ 2 ,/ 3 ) = F'(T,p',I 1 ,/ 2 ), similar to that for the equilibrium part 
the formulae for the complete stress tensor in the irreversible case: 

F° , results in 

g = -p T S + g- p T =Po+Pi> £“£ 0 + =r 

(23) 

p 0 = p 2 dF° /dp, p, = pp‘dF x / dp‘ ; 

(24) 

& 0 =2p[Ff(B-\/3lfS)-F 2 °(B~' -U3ifS)] . 

(25) 

d = 2 p[Ff (| - 1 / 37, S) - Ff (P - 1 / 3 1 2 £)] . 

(26) 



Here p T is the thermodynamic pressure due to compressibility of the viscoelastic solid, 
with equilibrium p Q and non-equilibrium p l components; <r is the thermodynamic stress 
deviator consisting of two, equilibrium d_ Q and non-equilibrium <5^ , deviator components. 


5 SIMPLIFICATIONS AND MULTI-MODAL GENERALIZATIONS 

Although only the theory with instantaneous elastic reply is considered in Section 3 and 
below, it still requires more simplifications to be useful in practice. We now consider 
several modeling arguments and corresponding steps introducing more simplicity in the 
theory. 

(i) The first step directed toward simplifying CE’s has already been made in eq.(l), which 
in fact is decoupling the free energy into a sum of equilibrium and quasi-equilibrium 
parts. Although there are some fundamental questions regarding the validity of this 
decoupling, it has been successfully used in many non-equilibrium theories. 

(ii) In spirit of the point (i), it is also assumed that it is possible to decouple the 
irreversible effects from the equilibrium response. It practically means that the kinetic 
tensor M in eq.(16) depends on the temperature and non-equilibrium elastic tensor b. 

This immediately reduces the evolution equation for tensor Mo the form: 

b+<p(T,b) = 0. (10a) 


2 3 ♦ • 

This form has successfully been used in the thermodynamic theories ' for viscoelastic 
liquids. Here in (10a) <p(T,b) is an isotropic tensor function of tensor b , whose properties 

2 3 

and specifications have been established in Refs. ' 

(iii) Decoupling the bulk and shearing properties. This can be justified in the case of 
large strains with relatively very small bulk deformations, which is the case of rubber 
elasticity. 

(iiia) Decoupling of the bulk and shearing thermodynamic properties. For the free 
energy this means the additivity: 


F = F° (7; If ) + F' (7; I ) ) = F b °(T; p) + F° (T\ I? , /* ) + (p‘ / p 0 )F l h (T ; p‘ ) + F sh (7; /, , / 2 ). 

(27) 

The stress tensor due to eq.(27) is therefore represented as: 


a = g 0 (T,B) + a k (T,b k ): 


ip 0 (T,p)+p k (T,p i )]S+g(T,B)+g k (T,b k ) (28) 


(iiib) Decoupling of the bulk and shearing relaxation properties. Eq.(22) with the use 
of simplification in (ii) (see eq.(10a)) can be represented as: 


d In p' / dt +[a(p")/d l b (T)](p i ) 2 dFfr / dp' (T,p‘) = d In p/dt (29) 

0 2 
b-b e-e b+ l/2[/3(I l J 2 )/8 l (T)] [b -1/3(7, -I 2 )b-S) = 0. 


( 30 ) 


Equations (27)-(30) impose the full decoupling of bulk and shear viscoelastic evolution 
equations. Due to the Second Law, the functions of/?') and /?(/,, / 2 ) are positive. They 

A A 

go to the unity in the equilibrium limit, when p 1 — /, — > / 2 ~ ^ 3 , and represent the 

strain scaling factors for corresponding relaxation times, which are similar in concept to 
the “material clocks” of the K-BKZ theory. The dissipative term in eq.(29) is written in 
general form, whereas in eq.(30) it has been presented in the form proposed and tested for 

polymer fluids and a gum rubbers 3 . The temperature dependent parameters 6 l b ( T ) and 
0 { (T) are the bulk and shear relaxation times, which occur in linear viscoelasticity. 

(iv) Simplified models for the heat conductivity tensor rc in (13) can be suggested when 
the tensor depends only on the equilibrium strain B, i.e. k_- k{T ,B). This dependence 

predicts that after unloading the heat conductivity is isotropic. Therefore the crucial test 
is measuring heat conductivity in a retardation experiment (i.e. after unloading). 

In the multi-mode generalization of CE’s for nonlinear solid-like viscoelasticity, we 
adopt all of the modeling simplifications discussed above for a single relaxation mode. 
We additionally advance the assumption that both the bulk and shearing relaxation times 
are well separated: 


0 k CD » GO; 6t CO » &L CO- (3 j ) 

Inequalities (31) which hold for viscoelastic liquids 3 give the opportunity to consider 
overall nonlinear properties as the sum of nonlinear viscoelastic properties for 
independent relaxation modes, where the free energy and stress (both bulk and shear 
components) are the sums of these for nonlinear relaxation modes. Additionally, the 
evolution equations of type (29), (30) also hold. For cross-linked rubbers, the fore 
mentioned closed set of nonlinear constitutive equations of differential type has the nice 
feature: they can model the properties of rubber-like materials in both the flow and solid 
regions with changing parameters during cure reaction, when keeping the same 
mathematical structure of equations, except for the equilibrium part. 

Additionally, some scaling approaches, like well-known time-temperature scaling of 
WLF type, could apply to both shear and bulk relaxations. Furthermore, the deformation 
scaling 24 developed for viscoelastic liquids, if applicable, could significantly simplify 
numerical computations for viscoelastic solids, too. 

6 STABILITY CONSTRAINTS FOR CONSTITUTIVE EQUATIONS FOR 
VISCOELASTIC SOLIDS 

The formulation of field equations for viscoelastic solids should obey some stability 
constraints. These constraints take into account the most fundamental properties of CE’s, 
related to the Second Law of thermodynamics, propagating various fast waves, and the 
fundamental condition providing the continuity for a solution of a problem along the time 
axis. Unlike CE’s for viscoelastic liquids, in the case of viscoelastic solids, these 
constraints are related only to the Hadamard stability (or the strong ellipticity, or the 
strong hyperbolicity for wave propagating, or well posedness). 



In the case of finite elasticity, the general (necessary and sufficient) condition of 
strong ellipticity is: 


A = A ijmn {B)x iyj x m y n >0. (32) 

The rank of four tensor A ijnm is defined for isotropic solids as: 


3<J, d 

{ d F ] 

An d 

f df) 

A ,JB)= " = 

- dh mn dh mn 

l P *\\ 


P " M 

) 


Here h is the Hencky strain measure. Condition (32) of strong ellipticity means the 

positive definiteness of the tensor A ijmn considered in terms of a bi-quadratic form. It was 
reduced in both incompressible 25 and compressible 26,27 cases to a set of algebraic 
inequalities imposed on the first and second free energy derivatives taken with respect to 
three basic invariants I k . In the case of the linear stability analysis of disturbances 
imposed on arbitrary, nonlinear, and generally time dependent stress-strain deformation 
field, condition (32) of strong ellipticity coincides with the condition of strong 
hyperbolicity: 


P& 2 v, v, = A ijnm ( B)v { k j v m k n > 0 . (34) 

Here Q. = co-k v is the frequency oscillations in the disturbance wave with a Doppler’s 
shift of the basic velocity field v = v(x,t) being taken into account, where h and v ; are 

respectively the components of the wave vector and the vector amplitude of the velocity 
disturbances. In the limit co — > °°, k — > °°one can consider a monochromatic plane wave 
of disturbances as being restricted to a local analysis of stability. It has been proved (e.g. 
see 28 ) that when a system is described by the quasi-linear set of PDE, the linear stability 
analysis coincides with more general method of characteristics. Positive definiteness of 
the form (34) means the possibility to continue in the positive direction of time axis, at 
any time instant, any solution of the equations for finite elasticity. If this is impossible, 
the characteristics are turned back, and near a turning point, their derivatives approach to 
infinity. Thus, Hadamard instability is a typical ill posedness accompanying by the blow- 
up instability. In the general case of compressibility, all possible types of waves of 
disturbances are allowed, including longitudinal waves. However, when an elastic 
material is treated as incompressible, i.e. p- const, the longitudinal waves of 
disturbances are forbidden. Thus because of the incompressibility condition, the wave 
and the velocity amplitude vectors are orthogonal: kpj } = 0 . 

When analyzing static stability, where infinitesimal (symmetric) strain disturbances 
£ are imposed on a given finite elastic strain B q , so that J? = J? 0 +£, one obtains 

(instead of the strong ellipticity condition of eq.(32)), a thermodynamic stability 
condition: 


A , =A. ( B)e e >0 

T ijmn y—' ij mn 


(35) 


Although the fourth-rank tensor Ai jmn is once again defined by eq.(33), it is easy to see 
that the thermodynamic condition of stability (35) is only a necessary condition for 
strong ellipticity (32). 

The necessary and sufficient conditions for strong ellipticity are too complex and 
awkward to be applied to every newly proposed CE. Therefore, several attempts have 
been undertaken to find non-trivial, yet simple enough, sufficient conditions of strong 
ellipticity. 

Renardy 29 proved that the CE’s for incompressible elastic and viscoelastic solids are 
strongly elliptic, if the strain energy function F = / r (7',/ 1 ,/ 2 ), which here depends only 
on two basic invariants, 7y and h, can be represented as monotonically increasing and 
convex function of -yJT^ and -<Jl, . More complicated sufficient conditions for strong 
ellipticity in compressible case have been recently derived in paper . 

7 REMARKS OF SPECIFICATION OF FUNCTIONS AND PARAMETERS 

This Section briefly surveys the various literature specifications of mechanical properties 
for rubber-like materials. More detailed description will be given in presentation and 
elsewhere. 

It has been documented 3 ' 32 that the bulk properties of cross-linked rubbers are well 
described by the van der Waals equation, where stability constraints can also be imposed. 
Although the most popular for numerical applications is the Mooney-Rivlin potential, its 
descriptive ability of shearing properties for the rubber materials is rather limited 33 . 

To the authors’ knowledge, the widest range in description of various experimental 
data for the cross-linked rubbers has been achieved using the Blatz-Sharda-Tschoegl 

(BST) potential 34 , which operates with the generalized strain measure B" . However, their 

model does not lend itself to computational efficiency, since these problems can be 
effectively treated only in the principal axes of tensor B . The stability constraints for this 

potential are unknown, mostly because of the generalized strain measure involved. 

Another elastic potential that has been recently proposed and well tested for 
viscoelastic liquids 35 can also be used for crosslinked rubbers. This potential is a 
combination of the Mooney-Rivlin and Knowles 36 potentials, widely used for 
calculations of large deformations of elastic solids. Using the Renardy sufficient 
condition of strong ellipticity, some sufficient stability criteria for this potential might be 
able to be established. 

Several recent attempts were also made to propose expressions forelastic potentials 
for very large deformations of elastic solids (e.g., see 37 ' 38 ), related to the finite 
extensibility of elastomer chains. These approaches used the inverse Langevin function, 
which is too complicated for solving complex problems. The easiest semi-empirical way 
to include this effect into consideration was proposed by Gent 39 , who was unaware that 
the same approach has been proposed for viscoelastic liquids (see Ref. 5 ) utilizing 
Warren’s approximation of the inverse Langevin function (“Warren spring” 40 ). This type 
of potential, depending only on 7, , might describe well enough very large elastic 


deformations. These usually occur in uniaxial extension or in situations close to that, 
which can be met in complex deformations with high-level concentrations of stresses. 
However, at relatively low strains, predictions of this potential are in contradiction with 
experimental data 33 , which clearly show the dependence of elastic stress on the second 
invariant 7, . This means that in low and intermediate regions of elastic strains, the 
potential 35 might be more preferable. It is possible to combine the potentials proposed in 
papers 35-36 for describing the large range of elastic deformations presented in rubbers. 

To describe the non-equilibrium nonlinear viscoelastic properties of crosslinked 
rubbers one has to specify the elastic potentials for the tensors b k and scaling factors 

presented in the evolution equations. The elastic potential that describes bulk non- 
equilibrium properties could be taken, once again, in the van der Waals form ’ ", and a 
potential for shear properties could be taken in the form used in Ref. 35 for viscoelastic 
liquids. The function /? , which describes the dissipative properties in shear evolution 
equations (36), could also be extracted from Ref. 35 , where it has been applied to 
viscoelastic liquids. The function a(p[), which describes the dissipative properties in 
bulk evolution equations (35), has yet to be determined via experimental studies and/or 
theoretical modeling. 

Finally, the effects of stress induced anisotropy in thermal conductivity found in 
papers 42-43 for uniaxial extension, could be easily generalized. This will be demonstrated 
in the presentation. 
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